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Abstract 

The positive weight next-to-leading-order matching formaUsm (POWHEG) is apphed to 
Deep Inelastic Scattering (DIS) and the related Higgs boson production via vector-boson 
fusion process in the Herwig-|— |- Monte Carlo event generator. This scheme combines 
parton shower simulation and next-to-leading-order calculation in a consistent way which 
only produces positive weight events. The simulation contains a full implementation of the 
truncated shower required to correctly model soft emissions in an angular-ordered parton 
shower. 
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1 Introduction 

The Large Hadron Collider (LHC) at CERN is designed to elucidate the nature 
of electroweak symmetry breaking in the Standard Model P-S], and in particular 
discover the Higgs boson. Once the Higgs boson has been observed and its mass 
determined, it will be crucial to measure the way it couples to gauge bosons and 
fermions |5l|6]. The most promising processes in which these couplings of the Higgs 
boson can be measured are gluon-gluon and vector-boson fusion. The former consists 
of a gluon-gluon partonic collision which produces the Higgs boson via a virtual top 
quark loop [7]. It has the largest cross section for Higgs boson masses less than 
~ 1 TeV and will be important for the measurement of the Higgs boson coupling to 
the top quark. 

Higgs boson production via vector-boson fusion (VBF) is a process in which two 
incoming fermions each radiate a or boson which then combine to produce 
the Higgs boson. VBF is expected to play a fundamental role in the measurement 
of the Higgs boson couplings to gauge bosons and fermions, because it allows for 
independent observation in different channels: if — >■ rr ^i9j, H — )■ WW [TO|ITT]. 
if -> 77 and ii ^invisible [T3l[Tl]. 

In order to calculate the Higgs boson coupling constants with sufficient accu- 
racy, next-to-leading-order (NLO) QCD corrections for the VBF process must be 
included. These corrections have been known for some time [T5] and are relatively 
small with iT-factors around 1.05 to 1.1. At next-to- leading-order, the theoretical 
prediction of the Standard Model production cross sections have an error of less 
than 10%. This accuracy is sufficient to compare predictions with upcoming LHC 
measurements, which will be performed with a statistical accuracy on the product 
of the production cross section and decay branching ratio reaching 5 to 10% [51,16]. 
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The theoretical uncertainties for the VBF process therefore do not significantly com- 
promise the precision of the coupling constant measurements. This makes the VBF 
process more attractive than Higgs boson production via gluon fusion, which has 
a K-factor larger than 2 and for which the uncertainties remain between 10 — 20% 
even after the inclusion of next-to-next-to- leading-order corrections [T6H23]. Nev- 
ertheless, stringent cuts are necessary to distinguish the VBF Higgs boson signal 
from the backgrounds. In particular, a veto on additional activity in the events, the 
central-jet veto, is often imposed to reduce the backgrounds. 

In order to study the effects of these cuts we must rely on Monte Carlo event 
generators which combine, usually leading-order (LO), matrix elements with parton 
showers and hadronization models to provide a fully exclusive simulation of an event. 
Traditionally leading-order matrix elements have been used in these simulations 
together with the parton shower approximation which simulates soft and collinear 
emission. However, in recent years a number of different approaches have been 
developed to improve the simulation of high transverse momentum, p^, radiatioiiEI. 

A number of approaches has been developed to provide a description of hardest 
emission together with a cross section which is accurate to next-to-leading-order. 
In the approach of Frixione and Webber (MC@NLO) [55t l56]. the parton shower 
approximation is subtracted from the exact next-to-leading-order calculation. This 
was the first successful systematic scheme for matching next-to-leading-order calcu- 
lations and parton showers and has been applied to many different processes [57T - [6^ . 
However, this method has two drawbacks: it generates weights which are not posi- 
tive definite and is implemented in a way which is fundamentally dependent on the 
details of the parton shower algorithm. 

These problems have been addressed with a new matching algorithm proposed 
by Nason, POWHEG (POsitive Weight Hardest Emission Generator) [65|[66]. which 
achieves the same aims as MC@NLO but produces only positive weight events. Al- 
though it is independent of the generator with which it is implemented, it requires 
the shower to have a particular structure: a truncated shower simulating wide angle 
soft emission; followed by the highest transverse momentum (px) parton emission; 
followed again by a vetoed shower simulating softer radiation. The highest pt emis- 
sion is generated separately using a Sudakov form factor containing the real emission 
piece of the differential cross section. The truncated shower produces radiation at a 
higher scale (in the evolution variable of the parton shower), and the vetoed shower 
at a lower scale, than the one at which the hardest emission is generated. The 
POWHEG method has been applied to wide range of processes [671 - 186 j^ . 

For the VBF process, as the central-jet veto is sensitive to additional radiation 
in the event, it is important to have an accurate simulation which gives both the 
NLO cross section for the process and provides an accurate simulation of additional 
QCD radiation. In this work we will describe the simulation of this process in 
Herwig++ [l3l|90] using the POWHEG approach. 



^ See Ref. [3S] for a recent review of the older techniques PS1H5] and techniques for improving 
the simulation of multiple hard QCD radiation [HHS^ . 

^ There has also been some work combining either many NLO matrix elements ^87^ or the NLO 
matrix elements with subsequent emissions matched to leading-order matrix elements |88[I89) with 
the parton shower. 
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As we develop new simulations it is important to validate them using experimen- 
tal data wherever possible. While obviously there is no data on the vector boson 
fusion process, Deep Inelastic Scattering (DIS) has many of the same features, in 
particular in both processes it is important that the parton shower algorithm pre- 
serves the momentum of the space-like vector bosons. In addition the wealth of 
data from HERA makes deep inelastic scattering important for the tuning of the 
phenomenological parameters in Monte Carlo event generators. There are existing 
simulations of the VBF [75] and DIS processes [76] in the POWHEG approach. 
However, due to different treatment of this class of processes in the angular-ordered 
parton shower in Herwig++ (the different kinematical reconstruction of these pro- 
cesses, preserving the virtuality of the t-channel gauge bosons, and generation of the 
truncated shower) and the experimental importance of this processes it is important 
to have a range of different simulations with different approaches. In addition our 
factorized approach makes the extension to other colour- singlet vector-boson fusion 
processes simple, requiring only a calculation of the leading-order matrix element. 

The rest of the paper is organized as follows. In Sect. |2]we recap the POWHEG 
formulae of main interest for our description. The calculation of the leading-order 
kinematics with NLO accuracy in the POWHEG approach will be discussed in 
Sect. El A brief description of the generation of the hardest emission within Her- 
wig++ will be outhned in Sect. [Hand in Sect. |5]we give details of the implementation 
of truncated and vetoed showers in the program. Our results will be described in 
Sect. [61 We finally present our conclusions in Sect. [T] 



2 The POWHEG method 



In this section we introduce the details of the POWHEG algorithm. The NLO 
differential cross section for a given N-body process, can be written within the 
POWHEG approach as 



dor = B{^B)d^B 

where B{^b) is defined as 



A«(0) + ^^§^;^An{kri<i>B, $/?))d<f r 
n['^B> 



J2Dj{<^b 



d$ 



B{^b) is the leading-order contribution, $b the N-body phase-space variables of 
the LO process, V{^b) a finite contribution including unresolvable, real emission 
and virtual loop pieces, $k are the radiative variables describing the phase space 
for the emission of an extra parton, R{^b, ^r) is the matrix element including the 
radiation of an additional parton multiplied by the relevant parton flux factors and 
Di{^b,^r) are the counter terms regulating the singularities in R{^b,^r)- The 
modified Sudakov form factor is defined in terms of -R($_b, as 



exp 



cm^^§^^Oikri<^B,<^R) 



Pt) 



(3) 
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3 Calculation of B{^b) 



where /cT($B,$i?,) is equal to the transverse momentum of the emitted parton in 
the soft and colhnear hmits. 

The POWHEG formahsm requires that the N-body configuration is generated 
according to B{^b)- The hardest emission in the event is then generated using 
the Sudakov form factor given in Eqn.[3l As is simply the next-to-leading- 

order differential cross section integrated over the radiative variables, it is naturally 
positive, and therefore leads to the absence of events with negative weights. 

If the parton shower simulation is ordered in transverse momentum we can sim- 
ulate the process by first generating the hardest emission and then evolving the 
process using the parton shower from the + 1 parton final state forbidding any 
emissions with transverse momentum above that of the hardest one. However, for 
shower algorithms which are ordered in other variables, for example angular ordering 
in Herwig++, the hardest pt emission is not generated as the first emission in the 
parton shower. So the shower must be reorganized into a truncated shower which 
describes soft emissions, at higher evolution scales than the highest pt emission, to- 
gether with vetoed showers which describe emissions at lower evolution scales which 
are constrained to be softer than the hardest emission [651166] . 

We implement the POWHEG algorithm according to the following procedure: 

• generate an event according to Eqn.[Tl 

• directly hadronize the small fraction of non-radiative events; 

• map the radiative variables parameterizing the emission into the evolution 
scale, momentum fraction and azimuthal angle, {qh, Zh,(j)h), from which the 
parton shower would reconstruct identical momenta; 

• consider the initial A^-body configuration generated from B{^b) and evolve 
the parton emitting the extra radiation from the default starting scale down 
to Qh using the truncated shower; 

• insert a branching with parameters {qh, Zh, 4>h) into the shower when the evo- 
lution scale reaches qh, 

• generate px vetoed showers from all external legs. 

This simple approach allows us to correctly generate wide-angle soft radiation using 
the truncated shower which was absent in some earlier simulations. 

3 Calculation of B{^b) 

At leading order DIS is described by the Feynman diagram shown in Fig.[T^, together 
with appropriate crossings of the quark line. The leading-order diagram for the VBF 
process is shown in Fig. [2^, together with appropriate crossings of the quark lines. 
In principle other contributions to the VBF process should be considered: diagrams 
with the exchange of identical outgoing quarks, the quark annihilation processes 
gg — )■ Z* — i- ZH and qq — )■ W^* — )■ W^H with hadronic decays of the vector 
bosons. However, colour suppression and the large momentum transfer in the weak- 
boson propagators make the contribution coming from these additional processes 



5 



9 




Fig. 1: Feynman diagrams contributing to deep inelastic scattering at 0{as)'- leading 
order (a), virtual (b) and real emission (c-f) corrections. 



negligible in the phase-space regions where VBF can be observed experimentally, 
i.e. with widely separated quark jets of very large invariant mass pi] . 

At 0{as), the contributions coming from amplitudes in which the gluon is at- 
tached to both upper and lower quark lines in the VBF process vanish because the 
weak boson has no colour charge. The Feynman graphs contributing are therefore 
only the ones shown in Figs.[2]D-[2f where for simplicity, we just show radiation from 
the upper quark line. 

The corrections to the DIS and VBF processes are therefore the same provided 
that we take into account the corrections to both of the quark lines in the VBF 
process. In this section, we show the analytical contributions to the next-to-leading- 
order differential cross section for DIS and VBF. Collecting the real emission cross 
section, described in Sect. 13. virtual and collinear contributions, briefly discussed 
in Sect. 13. 2[ can be calculated. We then discuss how it is sampled within 

Herwig++ in Sect. 13.31 
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3 Calculation of B{^b) 
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Fig. 2: Feynman diagrams contributing to Higgs boson production via weak-boson 
fusion at 0{as)'- leading order (a), virtual (b) and real emission (c-f) cor- 
rections. As matter of simplicity, we just show radiative corrections to the 
upper line in the gluon emission (c-d) and gluon initiated processes (e-f). In 
this case V = or Z". 

3.1 Real emission contribution 

The tree-level corrections to e'^e~ annihilation to hadrons can be written in a form 
in which QCD and electroweak parts exactly factorize [92]. This method was later 
generalized to any process in which the lowest-order diagrams contain a single quark 
line attached to a single electroweak gauge boson [37j. We adopt this approach for 
the calculation of the real corrections to DIS and VBF, based on the calculation of 
the correction to DIS in Refs. [33|[37]. 

Consider the QCD Compton process, shown in Fig.|3l where a quark q with 
momentum pi and a fraction Xb of the incoming hadron momentum, interacts with 
a current and boson-parton coupling u^, and scatters into an outgoing quark q' 
with momentum p2 and a gluon with momentum 



QCDC : + X ^ q'ip2) + gips) + X' . 



(4) 



3.1 Real emission contribution 
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Fig. 3: QCD Compton scattering, in which a quark interacts with an arbitrary cur- 
rent via boson-quark couphng u^. 



It is simplest to work in the Breit frame, in which the incoming parton for the 
leading-order process has four-momentum qi = ^(1; 0,0,1), the exchanged boson 
has four-momentum q = (0; 0, 0, —Q) and the scattered quark has four-momentum 
Q2 = ^(1; 0, 0, —1). In this frame the four-momenta of the real emission process are 

Pi = ^{-xi;0,0,-xi); (5a) 
P2 = ^i\f4 + x]_; x±cos4i, x_Lsin(/), -xs); (5b) 



Ps = y + -x± cos0, -x± sin 0, -X3); (5c) 

where the transfered momentum q = (0; 0, 0, —Q) = P2 + P3 — Pi and 

Xi = . (6) 

q ■ q 

Momentum conservation requires that X3 = 2 + Xi — X2 and 

™2 _ 1-^3 -^1 -^2) _ 2 (ry^ 
Xj_- ^_^2 •^2- {') 

In terms of these variables the cross section for the real emission process is 

_ 1 d(l) dXidX2-XiXBfi-XiXB,Q'^) Q'^IMqcdc? .^x 
''''^"°-4(2vr)2 2vr -xf Xsfixn^Q') \M^oiquq2W 

where is the momentum fraction of the quark in the leading-order process, 
lA^QCDCP and |A^lo('?i, 92)^ are the spin and colour averaged matrix elements 
squared for the real emission and leading-order processes respectively. 

Using the gauge introduced by the CALKUL collaboration [93j, the matrix ele- 
ment for the real emission process is 



I-MqcdcP = -jrT^^^r^^2 (^1 + (^2 + 4)^2) \M^o{qi, g2)r, (9) 

where 

p _ \MLo{n,r2)\^ 

xl + xl\M^o{qi,q2)r ^ ^ 
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3 Calculation of B{^b) 



with r2 = — and ri = r2 — q. 



The integration of the phase space is simpler if we use the variables: 

1 , , 

Xi = ; (11a) 

Xp 

x, = l- (lib) 

Xp 

so that the phase-space limits become xb < Xp < 1 and < < 1. In terms of 
these variables 

^2 _ 4(1 — Xp){\ — Zp)zp ^^2^ 

Xp 

The cross section for the real emission process becomes 

asCp d^ l^fi^^Q^) dxpdzp ^ 2/ 2 , 2nd^j 
= ^2^ xsfixB,Q^) il-xp)il-zp) + "-("^ + ^^'^ 

This allows us to treat the QCD Compton process as a correction to the leading- 
order quark scattering process. 
The boson-gluon fusion process, 

BGF : g{p^) + X ^ q'{p,) + q{p,) + X', (14) 

can be treated in a similar way. In this case the spin and colour averaged matrix 
element squared is 

\Mbgf\' = ^"1^."^^ ((^2 + 4)^2 + ixl + xl)Rs) , (15) 

[i — X2)[i. — XsjL^ 

where 

xj |.MLo(r3,r3 + g)p 
xi + xi |A^Lo(gi,g2)|^ ' ^ ' 

with rg = -P3/X3. 

Using the same change of variables as before the differential cross section is 

xb -f I xb_ /02 j 

Tnas dct) dXpdzp V ' ^ / / 2 / 2 , 2 x p . 2/2 , ^2NpN 

- -^2^zp{l-zp) xnf{xB,Q') ^ + ' + ^ + ■ 

(17) 

So far, our result gives the calculation of the BGF cross section, without any dis- 
tinction between quark and antiquark scattering. If we want to view Eqn.[T7] as a 
correction to a given lowest-order process, partons and antipartons should be treated 
equivalently. As the Zp = 1 singularity is associated with configurations that become 
coUinear to the lowest-order quark scattering process, while the Zp = singularity 
is associated with the antiquark scattering process, we can separate 

111 

r = - + z 18 



3.1 Real emission contribution 
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and rewrite the cross section as 

TrOs d(j) dXpdZp l^fi^) , 2/2, 2nd I 2/2, 2^D^ n n\ 

271 2n (1 - Zp) xbj[xb) 

with the corresponding term giving a correction to the antiquark scattering pro- 
cess [331I37]. 

Using these results we can rewrite the real emission corrections as 

daR^R{^B,^R)d^Bd^R, (20) 

where 

1 1 

with I e {QCDC, BGF} and 



Cqcdc = Cf, (22a) 
Cbgf = Tn, (22b) 

- p!|?!(i-^.Ki-^) ^ ^^^^^ + ^^^^^^^ ' ^^^^^ 

^BGF = ^i'l'^^Zln ^ ^ + 4)^2 + x^(x^ + xi)i?3) . (22d) 

XBjq[XB,Q^) (1 - ^p) 

The radiative phase-space element is 

d$i? = ^dxpdzpdcj). (23) 
The singularities in -R($s, $_r) are cancelled by subtracting 

D. = '^J^V,, (24) 
where Vi are the Catani-Seymour dipoles [91] : 



EBf(EB n2\ 2 , ^2 



^^QCDc - ^~rr^~~F;2Y7i — ^ttt^ — ^'^^^^ 

XBjqiXB, Q^) (1 - Xp)[l - Zp) 

|^/,(f^,Q')x2+(l-Xp)2 ^ ^ 

XBfgixB, Q^) Zp 

The contribution of the real emission processes to B is therefore 

B{^B)d^B ^^^^{Ai-Vr)d<^B,. (26) 

/e{QCDC,BGF} 

While these expressions are sufficient for the calculation of the real correction, 
in the case of DIS the leading-order matrix elements are simple enough that i?2,3 



10 



3 Calculation of B{^b) 



can be calculated analytically and integrated over the azimuthal angle, (f) and Zp to 
simplify the numerical integration of B. 
For DIS 



_ cos2 02 + A COS 62 {i - y/W^ sin 62 cos (p) + {i - \fW^ sin Q2 cos 

(27) 

where cos 6*2 = , ^ , sin^^a = , ^-^ ^ , £ = ^ — 1, vb = and is the four- 

momentum of the incoming lepton. A is related to the couplings of the fermions to 
the exchanged vector bosons. For the charged current process A = 1^ whereas for 
the neutral current process 

with r = (^^i^ and 

Cy, = ■ , ^ , f ^ - Q^ sin^ ^^^'j , (29a) 
sm fc'vi/ cos b'vi/ V ^ / 

Ca, = ■ . ^ . ^, (29b) 
sm c^vi/ cos o'pv' 2 

where m^: is the boson mass, 6w the Weinberg angle, Qi the fermion charge and 
Ts^i its weak isospin. 

The expression for can be obtained from that for R2 with the exchange 
A — > —A, 02 O3 and — )■ vr — 0. 

In this case the contribution to B is 

Bi^B)d^BdXp Yl ^^i^ (30) 

/e{QCDC,BGF} 

where 

ff/</(?'^')2 + 2£2-Xp + 3x/2 + ^£(l + 2xp) , 

^ Q') 1 + £^ + 2 (1 - 3i^) Xpjl - Xp) + 2Ai (1 - 2xp(l - Xp)) 

''''^ XBf,{xB,Q^) 1 + Ai + e 

(31b) 

3.2 Virtual contribution and collinear remainders 

The finite piece of the virtual correction is given by [91] 



day = ^^^VixB)Bi^B), (32) 

where the finite contribution of the 1(e) operator of Ref. [01] and the virtual cor- 
rection is 



VixB) = -^-^ + ^ln ^ + 2 ln(l - xb) In ^ + ln2(l - xb). (33) 



3.3 Sampling B within Herwig++ 
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The collinear remainders are 



ClCXcoll = 7^ —-D^'Pb) 



with the modified PDplfl 



2-n f{xB,fJ-F) 
1 dx. 



(34) 



Xb 



fg — ^f^F] A{x 



x< 



+ 



Xb 



Xb 



fq - XpfqiXB^jJ'F) 



BiXr, 



fiF C{X 



PJ I 5 



where fq and fg are the quark and gluon PDFs respectively, and 



T 

M^p) = 7^ K + (1 - In 



Q\l-x,) , ^Tf 



2 ~l~ 2^ X. 



pii 



BiXr, 



C{xp 



2 _^^Q'i^-Xp) 3 1 



X Xr. 



^4 

InXp — (1 + Xpjln — 



1 — X. 



up fl pXp 

The combined contribution of the finite virtual term and collinear remnants is 



da- 



2tt 



-V($b)5(<I>b) 



where 



with VixB,f^F) - ^"^"^'^^^ 



V{^B) = V{xB) + V{xB,flF), 

I 

f{xB,tJ.F) 



(35) 

(36) 
(37) 
(38) 

(39) 
(40) 



3.3 Sampling B within Herwig++ 

Using the results in the previous sections 



5f$ 



B 



5f$ 



B 



1 + 



2tt 



-V($b) 



+ 



2n 



[Aii<^B,^R) -'Dji^B,^R)]d^R 



■ (41) 



/6{QCDC,BGF} 

For convenience, the radiative variables {xp, Zp, 0} are transformed into a new 
set defined on the interval [0, 1], such that the radiative phase space is 

a unit cube. The variable Xp IS redefined as 



Xp = 1 — p 



1 

l-n 



(42) 



^ We write the modified PDF for a quark q, but a similar expression is valid for an incoming 
antiquark q. 
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4 The generation of the hardest emission 



where n is fixed, and p is the new variable with phase-space hmits 

0<p< (43) 

This change of variable has been made in order to guarantee numerical stability in 
calculating the integral of 1/(1 — Xp). A further transformation is needed to achieve 

P=(1-xb)'-%. (44) 

Finally, the variable (j) is easily redefined as = 

The sampling of B{^b) proceeds in the following way: 

1. generate a leading-order configuration using the standard Herwig+-)- leading- 
order matrix element generator, providing the Born variables $b with an as- 
sociated weight B{^b)', 

2. generate radiative variables, by sampling i?($B),which is parameterized 
in terms of the unit cube (xp, Zp, 0), and using the Auto-Compensating Divide- 
and-Conquer (ACDC) phase-space generator [95] : 

3. accepted the leading-order configuration with a probability proportional to the 
integrand of Eqn. ST] evaluated at {$_b, 

In order to treat radiation from both quark lines in the VBF process we randomly 

select one line which emits the radiation and generate events in $ jxp, Zp, 0|. The 

symmetry of the process then ensures that the correct statistical result is obtained 
by multiplying the correction term in Eqn.UT]by two. 



4 The generation of the hardest emission 

The hardest emission is generated using the modified Sudakov form factor, given 
by the product of A/j(pt) for each channel contributing; this is done replacing the 
ratio R{^b,^r)/B{^b) in Eqn.[3]with 

Moreover, we prefer to generate the hardest emission in terms of radiative variables 
$^(xx, Zp, (p) so that the 6'-function in Eqn. [3] simply gives x± as the upper limit of 
the integral and the modified Sudakov form factor, for the channel /, becomes 

~ / —;^d(j)dzp^-^8zp{l - Zp){l - XpfAi ) , (46) 

where yx™^^ is the maximum value for the transverse momentum. 

The radiative variables are generated using the veto algorithm, described in [96] . 
We use the upper bounding function 

91 = (47) 

for the integrand which is chosen so that gi can be easily integrated in {xj,, x™^^}. 
The generation procedure then proceeds as follows: 
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1. a;_L is set to x^^^; 

2. a new {xp, Zp, 0) is randomly generated according to 



A?e7(a;x) = exp ( - / " ^d^dz^a, ) , (4^ 



giving^' 



4 = — TT-^, (49) 

Zp = 7^2, (50) 

= 7^3. (51) 

3. if Xp < or Xp > 1, the configuration generated is outside the phase-space 
boundaries, set x^^^ to x± and return to step 1; 



4. if 

^ ^'""'-izpil - zp){l - XpfAj > 7^, (52) 



9i 27r 

the configuration is accepted, otherwise set xf^^ to x± and return to step 1. 



5 Truncated and vetoed parton showers 

The Herwig++ shower algorithm, [901 197], starts at an initial scale, given by the 
colour structure of the hard scattering process, and evolves down in the evolution 
variable related to the angular separation of parton branching products, q. The 
evolution is generated by the emission of partons in 1 — )■ 2 branching processes 
and each branching is described by a scale, q, a light-cone momentum fraction, z, 
and an azimuthal angle, 0. The latter parameters are used to uniquely define the 
momenta of all particles radiated in a shower. However, the Herwig++ approach 
generally requires some reshuffling of these momenta after the generation of the 
parton showers to ensure global energy-momentum conservation. 

{N + l)-body final states associated to the generation of the hardest emission 
are first interpreted as a standard Herwig-|--|- emission, from the A^-body configu- 
ration, specified by the branching variables {qu-, Zh-, (ph) ■ The POWHEG emission is 
performed as a single Herwig+-|- shower as follows: 

1. the truncated shower evolves from the default starting scale to q^, such that 
any further emission conserves the flavour of the emitting parton and has 
transverse momentum lower than that of the hardest emission; 

2. the hardest emission is forced with shower variables {(jh, Zh, 4>h)', 

3. the vetoed shower evolves down to the hadronization scale, vetoing any emis- 
sion with transverse momentum higher than that of the hardest emission. 



^ Here TZi defines a random number in [0,1 
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5 Truncated and vetoed parton showers 



The key feature of this approach is the abihty to interpret the hard emission 
in terms of the shower variables. In order to do this we first need to consider the 
treatment of processes with an initial-state final-state colour connection, such as 
DIS or VBF, in the Herwig++ parton shower. In these processes the momenta of 
the incoming and outgoing colour connected partons after the parton shower are 
first reconstructed from the shower variables as described in Ref. [90]. These off- 
shell momenta are such that energy and momentum is not conserved, so boosts 
are applied to the incoming and outgoing momenta such that the momentum of the 
virtual boson is preserved by the showering process. In DIS and VBF type processes 
the reconstructed momenta are boosted to the Breit-frame of the system before the 
radiation. We take pb to be the momentum of the original incoming parton and Pc 
to be the momentum of the original outgoing parton and Pa = Pc~ Pb, therefore in 
the Breit-frame 

p„ = Q(l; 0,0,-1). (53) 
We can then construct a set of basis vectors, 

=g(l; 0,0,1), n2 = Q(l; 0,0.-1). (54) 

The momenta of the off-shell incoming parton can be decomposed as 

gin = «inni + Ann2 + q±, (55) 

where am = An = and q_L = qin - ai^ni - /3in?^2- In order to reconstruct 

the final-state momentum we first apply a rotation so that the momentum of the 
outgoing jet is 

gout = ttout'^i + /3out'^2 + q±, (56) 

where /3out = 1 and the requirement that the virtual mass is preserved gives Oout = 
'^2n^nt ' ^^^^ I1. ~ ~p\- The momenta of the jets are rescaled such that 

/ , . An, out . fcn\ 

?in,out = "in.outfi^in.out'T'l + T "-2 + g±, (57) 

""in, out 

which ensures the virtual mass of the partons is preserved. The requirement that 
the momentum of the system is conserved, i.e. 

Pa = gL-gi'n = Q(o,o,-i;0), (58) 

gives 

Oin/^in - ftout^out = \, (59a) 



2' 

An _ /3out _ 1 
^in ^out 2 



(59b) 



Once the rescalings have been determined the jets are transformed using a boost 
such that 

boost / lp.r\\ 
gin,out > gin.ouf 



15 



In order to interpret the hard emission in terms of the shower variables we first 
calculate the momentum of the off-shell incoming, q'^'^'" , or outgoing parton, q'^'^'" , 
depending on whether we are dealing with initial- or final-state radiation. We then 
compute the boost into the Breit-frame of this system and construct the basis vectors 
ni^2 as before which allows us to determine the transverse momentum, g^, of the off- 
shell incoming parton. In this frame the momenta of the partons before the shower 
would be: 

P6 = |(1 + c; 0, 0, 1 + c); p,= |(1 + c; 0, 0, -(1 - c)). (61) 

The momenta of the off-shell partons before the boost required to conserve energy 
and momentum are 

qr = c^rn, + Ar^2 + g±, qr = <uTni + /3°rn2, (62a) 

where 

«r = Ar = J':'L ^ (63a) 

in 



rii ■ n2 2ni ■ n2(y 



-,/2 



Q-'iew anew Pc ' ^1 /^ou\ 

The inverse of the boost, which would be applied in the shower to ensure energy- 
momentum conservation, can then be determined and applied to all the incoming 
and outgoing partons. These momenta can then be decomposed in terms of the 
Sudakov basis used in Herwig++, allowing the shower variables {qt, Zh, (ph) to be 
determined. 



6 Results 

6.1 Deep Inelastic Scattering 

In order to test our implementation of the POWHEG approach for deep inelastic 
scattering we first compared the results from Herwig++ and DISENT |100j for the 
reduced cross section 

. ^ xQ^ dV 

27ra2y^ dxdQ2 ' ^ > 

where y = ^ and Y± = 1 ± (1 — y)^ and a is the fine-structure constant. The 
difference between the Herwig+-|- and DISENT results divided by the sum of the 
results is shown as a dashed line in the lower panels in Figs. [Hand [5] and is always less 
than one per mille. In addition Fig.H] shows the comparison of the Herwig++ result 
with the results from Ref. [HH] and Fig.E] shows the comparison of the Herwig+-t- 
result with the results from Refs. pS] and [99]. The excellent agreement with 
DISENT and the experimental data demonstrates that the generation of the Born 
variables and calculation of B is correct. In both cases the PDFs from Ref. |101] 
were used. 

In order to study the real emission we compare the results of Herwig++ with 
the measurements of the transverse energy flow in DIS from Ref. jl02j which are 
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Fig. 4: The e^p reduced cross section cr'^^^, as a function of x at fixed between 
200 GeV^ and 1500 GeV^ The experimental results of Ref. |98] are shown as 
crosses. The lower frame shows (Data — Theory) /Data and the yellow band 
gives the one sigma variation. The solid (black) line shows the Herwig++ 
result. 
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Fig. 5: The e^p reduced cross section a^^^, as a function of x at fixed between 
200 GeV^ and 1500 GeVl The experimental results of Ref. [M] for cj'^ ^ are 
shown as crosses and the results of Ref. [99] for cr^ ''as diamonds. The lower 
frame shows (Data — Theory)/Data and the inner (yellow) band gives the 
one sigma variation for a*^^^ and the outer (magenta) band the one sigma 
variation for 5"^ ^. The solid (black) line shows the Herwig++ result for a'^^^ 
and the dashed (red) line shower the Herwig++ result for ^. 
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Fig. 6: The inclusive transverse energy flow j^dE^/drj* at different values of x and 
for the low sample from |102] . The lower frame shows (Data — 
Theory) /Data and the yellow band gives the one sigma variation. 
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Fig. 7: The inclusive transverse energy flow j^dE^/drj* at different values of x and 
for the high sample from [T02j . The lower frame shows (Data — 
Theory) /Data and the yellow band gives the one sigma variation. 
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6 Results 



sensitive to the treatment of hard radiation in angular-ordered parton showers [33] . 
The comparison of Herwig++ with the low and high Q"^ samples from Ref. |102] are 
shown in Figs. E] and [71 respectively. In addition to the Herwig-|--|- result, with and 
without the POWHEG correction, we have included the result of the FORTRAN 
HERWIG [301I3T] and Herwig++ with a matrix element correction based on the 
approach of Ref. [33]. These results clearly show that without a correction to 
describe hard QCD radiation there is a deficit of emissions between 1 < r/* < 3 which 
is remedied by using either the POWHEG approach or a traditional matrix element 
correction. In general the POWHEG approach gives slightly less radiation than 
the matrix element due to the Sudakov suppression of radiation which is neglected 
in the matrix element correction approach and is in the best agreement with the 
experimental results. In these plots we have tuned the mass parameter for the 
splitting of soft beam remnant clusters in Herwig+-|- to 0.5 GeV from the HERWIG 
value of 1 GeV. The transverse energy flow in DIS is most sensitive to this parameter 
and the original HERWIG value was tuned to older transverse energy flow data. 

6.2 Higgs Boson Production via Vector Boson Fusion 

The typical feature of the VBF process at hadron colliders is the presence of two 
forward tagging jets. At leading order, they correspond to the two scattered quarks 
in the hard process and their observation, together with the properties of the Higgs 
boson decay products, is vital for the suppression of backgrounds [HHH] . The tagging 
jet distributions must be known precisely to gain a good estimate of the Higgs 
boson couplings: comparison of the Higgs boson production rate with the tagging 
jet cross section, within cuts, determines the Higgs boson couplings, [5l[6], and the 
uncertainties of the measured couplings are determined by the theoretical error of the 
cross section. At next-to-leading-order, the tagging jet distributions are sufficient 
to estimate size and uncertainties of the higher order QCD corrections, because 
the Higgs boson does not induce spin correlations in the phase space of its decay 
products. 

A detailed analysis of jet distributions has been realized in the present work and 
the results are shown in this section. A preliminary step has been the validation of 
the B function, by comparing the NLO differential cross section as function of the 
rapidity of the stable Higgs boson given by Herwig++ and VBFNLO, as shown in 
Fig.e 

However, at next-to-leading-order we can either encounter two jets, with one of 
them composed of two partons (recombination effects), or three jets corresponding 
to well-separated partons. As for LHC data, an algorithm to select two tagging jets 
is needed. There are two possibilities [91] : 

1. pT-method: the two tagging jets are the two highest pt jets in the event; 

2. E-method: the two tagging jets are the two highest energy jets in the event. 

We follow the pr-niethod and jets are defined according to the algorithm 
by using the Fast Jet package |lU3j . Cuts need to be chosen to reduce the effect of 
backgrounds and we follow the ones introduced in Tagging jets are required to 
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Fig. 8: Rapidity of the Higgs boson at NLO. Results from Herwig++ (solid black 
line) are compared to the one from VBFNLO (dashed red line). 



have transverse momentum, p^, and rapidity, yj, fulfilling the following cuts: 

Pt > 20 GeV, \yj\ < 4.5. (65) 
Moreover we generate the Higgs boson decay in r+r^ isotropically and require 



that the produced leptons have transverse momentum, px 
ity, r]^{+,-), so that 



A+.-)' 



PT^^^ _^ > 20 GeV, 



and pseudorapid- 



(66) 



In addition, we require that jet-lepton separation in the rapidity-azimuthal angle 
plane satisfies 

Ai?,.,(+,-) > 0.6, (67) 
and the tau leptons to lie between the two tagging jets in rapidity 

yj,mm < Vt(+-') < yj,mm- (68) 

Backgrounds to VBF are significantly suppressed if the two tagging jets are well 
separated in rapidity; therefore, we require 



(69) 



The factorization and the renormalization scale are chosen to be equal to the mass 
of the Higgs boson, tuh = 120 GeV. The other relevant electroweak parameters are 

Mw = 80.3980GeV, Mz = 91.1876GeV, a^^ = 0.007556, sin^^w = 0.222646, 

(70) 

and the weak coupling is computed as g = ■^4aem/sin6'w- The parton distribution 
functions are chosen to be the CTEQ6M set [10411105] . 
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6 Results 




Fig. 9: Difference of rapidity (left panel) and higher pt (right panel) distribution of 
the two tagging jets. In the left panel we have excluded the cut in Eqn. |69l 



The analysis provides the comparison of distributions for POWHEG implementa- 
tion (solid black curve) and LO simulation (dashed red curve) of Herwig++ parton 
shower together with VBFNLO NLO differential cross section (dotted blue curve). 

In Fig. [9] we present the differential cross section as a function of the rapidity 
separation and higher px of the two tagging jets. In the left panel we have excluded 
the cut in Eqn. |69l The cross sections show a peak at a rapidity of around 5 in the 
left panel and a transverse momentum of 70 GeV in the right panel. The POWHEG 
implementation leaves the cross section with respect to the px of the hardest jet 
unchanged, while it modifies the difference of rapidity distribution respect to the LO 
simulation of the Herwig++ parton shower: the peak is slightly lower and shifted to 
a higher value of the rapidity difference. 

In Fig. [To] we plot the cross section with respect to the transverse momentum 
(left panel) and rapidity (right panel) of the softer of the two tagging jets. The 
transverse momentum distribution shows a peak around 30 GeV and the rapidity 
around 2. The Herwig++ shower provides a similar description at LO and NLO 
accuracy. 

The transverse momentum and rapidity distributions of the third jet are plotted 
in Fig. [TT]in the left and right panel respectively. As would be expected here we see a 
harder spectrum for the third jet in the POWHEG approach which is now simulated 
using the real emission matrix element rather than the shower approximation. 

In these plots we see that the Herwig++ results lie below the fixed NLO results 
for distributions involving the two leading jets as a result of the subsequent parton 
shower, unlike the results of Ref. [75j where there is little difference between the 
POWHEG and fixed order results. This difference exists at both leading order and 
in the POWHEG approach in Herwig++ and is a result of the different shower algo- 
rithm and kinematic reconstruction in Herwig++. Given the excellent description of 
the related DIS data, compared with the previous HERWIG shower algorithm, this 
is an important difference in the two approaches and worthy of further study. 



6.2 Higgs Boson Production via Vector Boson Fusion 



23 



Herwig++ POWHRG 
Herwig++ T.O 
VBFNLO 




^ 0.0004 — 



Herwig + + POWHEG 

Herwig + + T.O 

■ VBFNT.O 




Fig. 10: Distributions of smaller transverse momentum, p^"^^^, (left panel) and 
smaller rapidity, yj^minj (right panel) of the two tagging jets. 



6.0-10~= — 



3. 0-10"'= — 



- Herwig++ POffHEG 
Herwig + + LO 
VBFNLO 



0,000125 — 





Hcr-w-ig + + POWHhIG 

--- Herwig + + LO 
■■ ■ VBFNLO 



Pt^' [GeV] 



Fig. 11: Transverse momentum distribution of the third jet (p^^) is plotted in the 
left panel and the rapidity distribution of the same jet {yj.^) in the right one. 
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7 Conclusion 

In this work the POWHEG NLO matching scheme has been implemented in the 
Herwig++ Monte Carlo event generator for DIS and Higgs boson production via 
vector boson fusion. For both hard subprocesses the B function has been calculated 
following the general approach provided in [36] which enables other colour-singlet 
VBF production processes to be easily included. The simulation contains a full 
treatment of the truncated shower, which is needed for the production of wide angle, 
soft radiation in angular-ordered parton showers. 

For the DIS implementation we find the cross sections to be in good agree- 
ment with the experimental results of Refs. [98ll99l[T02] . Our results show that the 
POWHEG approach correctly populates the so-called dead zone, as it appears in the 
transverse energy flow distributions at high Q^. 

For the VBF Higgs boson production we have shown different jet distributions 
after imposing typical cuts which are required to remove the effects of backgrounds. 
We find that the POWHEG implementation does improve the rapidity separation 
distribution for the two tagging jets and of pt and rapidity for the third hardest 
jet, while it mainly leaves the other distributions of the two tagging jets unchanged 
respect to the LO simulations within the Herwig++ parton shower. 

The lack of data prevents us from comparing the jet distributions with experi- 
mental results. However, the present work together with the Higgs production via 
gluon fusion and the Higgs-strahlung simulations, which were already implemented 
in Herwig++ 2.3 [73], will provide an important tool for analyzing the upcoming 
results at the LHC 
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